KSC VAB Aeroacoustic Hazard Assessment 


Justin M. Oliveira 1 * * , Sabrina Yedo : , and Michael D. Campbell 
NASA. Kennedy > Space Center , FL. 32899 

Joseph P. Atkinson 4 

ASRC Aerospace Corporation. Titusville, FL, 32780 

NASA Kennedy Space Center (KSC) carried out an analysis of the effects of aeroacoustics 
produced by stationary solid rocket motors in processing areas at KSC. In the current 
paper, attention is directed toward the acoustic effects of a motor burning w ithin the Vehicle 
Assembly Building (VAB). The analysis was carried out with support from ASRC Aerospace 
who modeled transmission effects into surrounding facilities. Calculations were done using 
semi-analytical models for both aeroacoustics and transmission. From the results it was 
concluded that acoustic hazards in proximity to the source of ignition and plume can be 
severe; acoustic hazards in the far-field are significantly lower. 


Nomenclature 

A area, nr 

A w A-weighted sound pressure level correction, dB 

a speed of sound, m/s 

b line support spacing, m 

C model coefficient 

c phase velocity’, m/s 

d diameter, m or interior wall spacing, m 

Dl directivity index, dB 

F thrust or force, N 

f frequency, Hz 

L length, m 

L w overall sound pressure level, dB 

M Mach number 

m mass per unit area, kg/nr 

m mass flow rate, kg/s 

OASPL overall sound pressure level, dB 

P, p pressure or overpressure. Pa 

R gas constant. J/kg-K or distance to receiver from plume source, m 

R a A-weighted frequency sensitivity weighting function, W 
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total absorption surface area, nr 

Strouhal Number 

sound pressure level, dB 

plume path length from nozzle to source, m 

temperature, K 

transmission loss, dB 

thickness, m 

velocity, m/s 

sound power, W or weight, lb 
width, m 

overall acoustic sound power, W 
plume core length, m 
absorption coefficient 

geometric angle from receiver to local plume slice flow direction, radians 
effective angle corrected for Strouhal Number effects, radians 
bandwidth of frequency band, Hz 
plume slice length, m 

acoustic efficiency, decimal or wall/window-panel material loss factor 

specific heat ratio 

wavelength, m 

Poisson's Ratio 

density, kg/nr 

modulus of elasticity, psi 

incidence angle, radians 

transmissivity 

Subscripts 

ambient 
point b 
nozzle/bore 
chamber or critical 
point d 

nozzle exit or point e 
frequency or point f 
incident 

outer sheet of sandwich 
plume slice 

plume slice, frequency 

reflected 

total 

nozzle/bore, blast-origin, or fundamental 
lower critical 
higher critical 
sonic condition 


I. Introduction 

r J^HE Vehicle Assembly Building (VAB) was built to allow processing of Apollo Program launch vehicles in the 
1960's and later was adapted to process the Space Shuttle launch vehicles. Since then the VAB has served 
effectively as the location to assemble and mate solid rocket boosters (SRB) to the vehicles External Tank (ET) and 
Orbiter. The VAB will again be adapted to process new human and non-human rated launch vehicles in support of 
NASA’s strategic plan for space exploration. 


Previously, for the Space Shuttle Program (SSP), the Kennedy Space Center (KSC) Safety Office determined a safe 
separation distance (SSD) of 1,315 feet for the inhabited building distance (1BD) [7]. The calculation used the 
Department of Defense (DoD) explosive safety standards (sometimes referred to as the weight based approach). The 
SSP safe separation distance for the VAB, derived using the DoD weight based approach, allows for processing of 
up to four complete Shuttle SRB's containing approximately 4.44 million pounds of solid propellant. The currently 
proposed Constellation Program (CxP) launch schedule could require assembly and storage of nearly 13 million 
pounds of hazard classification 1.3 propellants; up to eight 5-segment ARES-V boosters to be stored in the VAB 
simultaneously. Using the weight based approach for eight ARES-V boosters would result in a safe separation 
distance of 1,810 ft. This distance extends beyond a number of high occupancy facilities in the area around the 
VAB. The NESC concluded that the DoD weight based approach was overly conservative for SRMs with a 1.3 
hazard classification. NASA Engineering and Safety Center (NESC) recommended that an alternative approach 
based on calculating the threat due to each hazard component (heat flux, toxics, acoustics, etc.) be pursued [ 1 ]. 

In the event of an inadvertent ignition, several hazardous components could either prevent or restrict egress of 
personnel and potentially impact the surrounding facilities. A methodology development as well as an assessment of 
the noise levels in proximity of the VAB was required to support the alternative approach recommended by the 
NESC. Acoustic SSD's as well as levels of resistance nearby facilities have from acoustic energy was addressed. 
Several scenarios were considered and later in the study the case set was bounded using results from parallel studies. 
These parallel studies helped narrow the VAB building configuration and ignition propagation as a result of an 
inadvertent ignition. The emphasis of this study was placed on the far-field sound hazards to personnel and the 
effectiveness of surrounding facilities as safe heavens from acoustic energy. 

Prediction of aero-acoustic noise as well as transmission into nearby facilities was performed using semi-analytical 
models. The facilities modeled were those in immediate proximity to the VAB; the Launch Control Center (LCC) 
and Orbiter Processing Facility U 3 (OPF-3). An exposure criterion for personnel was developed to establish safe 
separation distances. All methodologies and findings were examined by a NASA peer review' committee. 


II. Modeling Assumptions 


Currently, CxP would reuse the VAB with minor modifications to support processing of the Ares-l/V vehicles. The 
VAB is the fourth largest building in the world by volume. It consists of 4 High Bays used to process NASA launch 
vehicles, each assembled on a Mobile Launch Platform (MLP). Several platform levels are used to perform 
processing operations and will too be reused and/or modified. Each HB is separated by concrete walls (lightly 
reinforced) and large steel columns which provide the main support to the structure. 



Figure 1: NASA Ares-l/V Flight Vehicles and the VAB 

The shell is made up of hundreds of aluminum “p unc h-ouf' panels which are designed to fail above approximately 
100 psf (meant to protect the facility from large pressure loading during severe weather). It was shown from CFD 
models that the pressure loading from ignition over-pressure (IOP) and plume impingement (resulting from an 
inadvertent ignition) would knock over internal walls and fail all VAB paneling. In summary, shortly after ignition 
only the VAB ceiling, main supports, and ML’s would remain intact; it was assumed the vehicles remain intact. This 
assertion was used in the acoustic modeling and helped to simplify the analysis such to make the acoustic fields 
diffuse into a fourth space. 


Though a fully processed vehicle would certainly become propulsive if an inadvertent ignition occurred, an un- 
capped booster would not. This is because the chamber pressure isn’t sufficient enough to create enough thrust. To 
simplify the analysis, it was assumed that all burning boosters were fixed in place in the VAB. Furthermore nozzle 
gimbling and motor pitching was ignored. 


Furthermore it was assumed that plumes emanating from boosters and impinging on the VAB floor/ceiling spread 
symmetrically. This assumption made the acoustic calculations easier because the fields could be treated as two- 
dimensional. Also, in scenarios were adjacent boosters (on the same MLP) burned together (and in the same 
direction), their plume acoustic powers were combined and an equivalent plume was computed. This simplification 
is typically done in aeroacoustics since the plumes become conjoined if in close proximity. 


Sound absorption (dissipation) from the atmosphere is dependent on temperature and relative humidity. The 


atmosphere at KSC is fairly humid year round (typically 40-95%) and temperature can vary between extremes 


(typically 35-105° F). Above approximately 4 kHz, temperature and relative humidity become competing effects. 


Atmospheric Absorption (Orlando Florida) 



Figure 2: Atmospheric Sound Absorption versus Relative Humidity (for various frequencies and temperatures) 


The dissipation from atmospheric absorption can be significant above 5 kHz. For this study the absorption was 


ignored to be conservative and simplify the analysis. Furthermore, ground reflections other than the initial plume 


deflection were ignored. 


III. Theoretical Background 

The computational modeling of aeroacoustics from plumes is difficult because of the small time steps and large flow 
times required to resolve spectra (10 kHz - 10' 4 seconds). Furthermore, models used to convert the unsteady- 


turbulent fields to acoustic sources are computationally expensive. Alternatively, semi-empirical methods have been 
used that combine theoretical equations with empirical correlations. The method used in this study breaks the plume 
into several slices which are converted to monopole sources (see figure below). 



Figure 3: Plume conversion into acoustic sources [2] 

The sources are then characterized and sound is radiated inside the domain. This method has proven to give good 
agreement with measured data. Most recent, comparisons with data from the Ares-IX launch showed phenomenal 
agreement (see a summary below). 


Proximity 

Spectrum Error 

OASPL Error 

Near-field; ML/Launch Pad (0-500') 
Mid-field; Launch Site (1000-1400') 

±0.5 dB 

+ 1 dB 

Far-field; adjacent pad/VAB (5000-10,000’) 

±5dB 


Figure 4: Model error determined from Ares-IX launch test data 


IV. Exposure Criteria 

Several sources stating the sound pressure levels associated with human exposure exist. OSH A standard 1019.95 [3] 
provides exposure limits for industrial areas based on exposure time. Because capped and uncapped boosters can 
bum for long periods of time (> 5 minutes), using the OSHA table would result in a 115 dBA exposure limit. The 
intent of the OSHA table was for use in industrial environments and not for inadvertent ignition-type events. 
Furthermore, use of the OSHA value would result in a large SSD (overly conservative). Hearing damage will largely 
affect a person's ability to egress because of the ear's sensitivity. 


A-Weighted Correction 



Figure 5: A - Weighted Correction [4] 

Ruptured ear drum(s) will impede a person's balance and ability to egress; this occurs around 160 dB [4]. The ear's 
threshold of pain occurs around 140 dB; intense nausea is typical at levels around and above 158 dB [4]. 

The personnel exposure limit used in this study was 140 dB (corresponding to the threshold of pain). At these levels 
personnel will experience a large degree of permanent and irreversible hearing damage but not ear drum rupture; 
thus allowing the ability to egress. Because of the log(R 2 ) nature of acoustic sound propagation, revising the 
criteria upward will drastically reduce the SSD. Also, OASPL levels are measured at a head level of 5' 6" and 
revision of this target height to a larger value would also drastically reduce the SSD. 

V. Calculation Methodology 

Acoustic predictions were performed using a modified Eldred's 2nd Method. The plume is broken into several slices 
and given a unique spectrum. The sound pressure from all the slices (at each octave band frequency) is summed to 
determine the net sound pressure at a specific point. This is done first by solving for the plume metrics using the 
chamber conditions. Using provided chamber conditions and nozzle geometry, nozzle exit conditions were found 
assuming ID compressible isentropic flow [6]. Assuming the flow to be choked at the nozzle throat, mass flow rate 


is found. 


m = 


Pc^_ 

jTc 


R c \y c + l) 


(y c +i)/(yc-i) 


Eq. 1 


The exit Mach number is found using the area-Mach number relation. 




(Vc+DAXc-i) 


Eq. 2 


The equation has two solutions; subsonic and supersonic condition. Since the chamber pressure is such that it’s 
above the critical pressure ratio, we assume the supersonic solution. 


Area-Mach Number Relation 



Figure 6: Area-Mach number relation (for various specific heat ratios ) 
The exit plane pressure is then found using the Mach number. 


= P e (l + \-Mj) 


_ i v -y c /(yc-i) 


2 ‘ 

Knowing the exit Mach number, the exit plane temperature is then found, 

T e = T c (\ + Y -^-Ml) 

With the exit temperature found, exit plane speed of sound is found. 


Eq. 3 


Eq. 4 


— V Yc^c^e 


Eq. 5 


The exit plane velocity is found then using the Mach number and speed of sound. 


Thrust is then found to be [9], 


F = mV e + ( P e - P amb )A e Eq. 7 

Plume core length is found using the exit Mach number-diameter correlation [2], 

^ = 3.45(1 + 0.38 M e ) 2 Eq. 8 

d e 

This correlation was derived from various experimental data sets of various plume types; it was a least-square-type 
fit to that data. 



Figure 7: Non-dimensional plume core length 


In order to calculate the overall sound power, an acoustic efficiency must be used. With direct floor impingement, 
the efficiency is found using the data in SP-8072 in conjunction with the nozzle geometry and placement relative to 
the VAB floor and/or ceiling. 
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Figure 8: Acoustic efficiency of various plume deflections [2] 


The overall sound power is then found using the efficiency factored against the mechanical power [4], 

W 0A = T]FV e Eq. 9 

The overall sound power can then be translated into an overall sound pressure level [2], 

L w = \0log{W OA ) -I- 120 Eq. 10 


Once the plume metrics are known, the plume is discretized into several slices and Strouhal numbers are found for 
each plume slice using the flow variables and choosing a frequency range (typically 1 0°- 1 0 4 Hz), 


S n (s) = 


fS0e 

^e^-amb 


Eq. II 


The normalize overall sound power of each slice is found, 


101og(x t W ' (s) / lVo J = 10 ^ 
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Eq. 12 


where Cj, c 2 , c 3 , c 4 , and c 5 are correlation coefficients (omitted for SBU considerations). The strongest source 
typically sits between 1.5-2 core lengths from the exit plane (shown in figure below). The core (inviscid) region 
terminates (x,) when the jet-edge shear layers collapse and form a large turbulent zone in the plume. This is also the 
cause of the log reduction in axial velocity along the plume axis. 



Figure 9: Normalized overall sound power along, plume 
The overall acoustic power for each slice is found [2], 


Lw.ps = 10 log 


x t W(s) 


+ L w + 




Eq. 13 


(m.'Vvm i^fi) = io/op ( , 


Eq. 14 


Knowing the Strouhal number range, the normalized power spectrum is found for each plume slice, 

r V e Gamb,\ «/%,__/ ^l^n (^) 2 

where C|, c 2 , c 3 , and c 4 are correlation coefficients (omitted for SBU considerations). An important caveat to these 
empirical correlations is that their coefficients are dependent on the propulsion type (chemical, nuclear, etc). 


Normalized Acoustic Power in Plume 
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Plume Slice Normalized Acoustic Power Spectrum 
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Figure 10: Normalized power spectrum of plume slice 
The power spectrum for each plume slice is then found [2], 

L w . ps , r = 10 log (w(f,s)/W(s){^22t]\ + L w , ps - 10 log + lOlog (A/) Eq. 15 

The directivity indices are found by determining the true angles the receiver location makes with each plume slice. 
The true angles are then converted to equivalent angles correcting for Strouhal number effects, 


Peff =P~ Ci log (c 2 S n (s)) Eq. 16 

where C| and c 2 are correlation coefficients (omitted for SBU considerations). Knowing the effective angle, 
directivity index is computed. 


Difarr) = io log 


Cl +c 2 cos 4 (/? ( , /7 ) 


- c 9 log(S„(s)) - c i0 


Eq. 17 


[(1 - c 3 cos(J} erf )) c < + c 5 ] r "[l + c 7 exp(c e p efr )\ 
where c h c 2 , c ? , c 4 , c 5 , c 6 , c 7 , c 8 , c 9 , and Cj 0 are correlation coefficients (omitted for SBU considerations). The 
average (over a typical Strouhal number range) directivity of chemical rockets it's approximately 40° (see figure 


below). 


Plume Source Directivity Index 
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Figure 1 1: Acoustic source directivity index (at various Strouhal numbers) 

Again, the coefficients vary for different propulsion types. The mean directivity for various propulsion sources is 
shown in the figure below. 



Figure 12: Directional characteristics for various propulsion types [2] 


The frequency dependant sound pressure level from each plume slice for a receiver is [2], 


Eq. 18 


SPL ps , f = L wps j - 10 log (R 2 ) - 11 + Dl{p eff ) 

A logarithmic summation over all plume slices of the above will yield the sound pressure spectrum for a receiver 


[ 2 ], 



Eq. 19 


The overall sound pressure level for a receiver can be found from another logarithmic summation over the entire 
spectrum [2], 


OASPL = 10 log 



Eq. 20 


Corrections to the sound pressure spectrum for ear sensitivity can be done by including the A-weighted term with 
the SPL psf equation. 


SPL f . Aw = SPLf + A w {f) 

The A-weighted correction for the human ear is frequency dependant [4], 

12200 2 / 4 


RM) = == 

( / 2 + 20.6 2 ) V (/ 2 + 107.7 2 )(/ 2 + 737.9 2 )(/ 2 + 12200 2 ) 

AM) = 2.0 + 20log(R A (/)) 

Summation of SPLu w will give the OASPL in dBA. 


Eq. 21 


Eq. 22 
Eq. 23 


Corrections to the sound pressure spectrum (normal not A-weighted) for diffraction with a wall can be done by 
including the diffraction term with the SPL ps r equation. An empirical correlation used to determine the diffracted 
gain in the acoustic spectrum (see figure below). 
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Figure 13: Diffracted gain versus non-dimensional characteristic length [5] 


The non-dimensional characteristic length can be converted to frequency by. 



The acoustic spectrums are then corrected for diffraction by applying the correlation. 


SPL f>diff = SPLf + DIFFtf) 


Eq. 24 


Eq. 25 


The diffracted acoustic spectrums are then applied to all walls of the facility. 

/1\ /i4(l-a)\ 

SPLf. inner = SPL fdiff - 10 log 10 (r) + 10 log, 0 ( — — — j Eq. 26 

In the above equation, the second term is the transmission loss (TL; from wall and windows) and third is the 
absorption loss (from internal carpeting, cubicles, etc). 


The average transmissivity is found by taking an area-weighted average for all walls. The transmission loss is 
depended on the wall fundamental and critical frequency. 



m(l2(l - v 2 )) 

J w 


Eq. 28 


The transmission loss is also largely dependent on the wall type (single or double panel wall). For a single panel 
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Eq. 29 
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It is important to note that the method conservatively assumes there is no transmission loss below the fundamental 
frequency of the panel. This was done primarily because specific treatment of this region was not taken into account 
being that the sound field is diffuse [11]. 
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Figure 14: Typical transmission loss profile for a single panel wall 
For double wall constructions, the analysis methodology is significantly different. First, the calculations for the 
method must be done in metric units. This is due to several simplifications that were included in the development of 


the procedure by Bies and Hansen [11]. Secondly, since double wall constructions can contain an air space, the 
equation for the fundamental frequency changes to account for coupling effects [8], 


/ o = 80 


(m, + m 2 ) 


Eq. 30 


d(m 1 m 2 ) 

Once this frequency is determined, it then becomes necessary to determine the critical frequency of both wall panels 

[ 8 ], 


f<=x.2 2n 


m,. 2 (l2(l -(v,. 2 ) 2 )) 


Eq. 31 


£1, 2(^.2) 

where if any wall panel is made of thin, equal-thickness sheets of material sandwiching a lightweight core, the 
critical frequency equation changes to [8], 
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Eq. 32 


It is important to note that the subscripts 1 and 2 in these equations correspond to the panel with the lowest critical 
frequency and highest critical frequency, respectively. The general transmission loss spectrum for a double panel 
construction is shown in the figure below. 


TL(ciB) 



Figure 15: Typical transmission loss profile for double panel wall 
For a double panel wall the transmission loss is found from. 
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Eq. 33 
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An important aspect of the double wall calculation method is that linear interpolation is required to solve for the 
transmission loss at the frequency values between the known points. This is due to the fact that there are several 
possible combinations of equations that need to be used, depending on the boundary conditions of the wall/window 
being analyzed. An example of this is when calculating the transmission loss at point B on the spectrum. If the space 
between the two panels comprising the double wall construction contains no sound absorbing material, the value of 
B| is used at point B; otherwise, the value at point B is the larger of B| and B 2 . 


The average absorption coefficient is found by taking an area-weighted average of ail individual absorption sources, 


— _ a i A i+ a i+i A i+i+'+ a n A n 
A i +A i+i + -+ A n 


Eq. 34 


Typical values for absorption coefficient are shown in the table below. 


ACOUSTIC ABSORPTION COEFFICIENTS 

MATERIALS 

125 HZ 

250 HZ 

500 HZ 

1000 HZ 

2000 HZ 

4000 HZ 

Painted Concrete: 

0.1 

0.05 

0.06 

0.07 

0.09 

0.08 

Carpet: 

0.01 

0.05 

0.10 

0.20 

0.45 

0.65 

Window Glass: 

0.18 

0.06 

0.04 

0.03 

0.02 

0.02 

Linoleum Floor: 

0.02 

0.03 

0.03 

0.03 

0.03 

0.02 

Ceiling Tiles: 

0.30 

0.51 

0.80 

0.85 

0.78 

0.66 

Gypsum Board 

0.29 

0.10 

0.06 

0.05 

0.04 

0.04 

Removable Walls 

0.14 

0.10 

0.06 

0.04 

0.04 

0.03 


Figure 16: Typical values of absorption coefficient [10] 


VI. Sample Scenarios 

Since CxP is the current NASA program, a sample case for the larger Ares-V heavy lift vehicle follows. In the 
scenario, an Ares-V vehicle is in the final stages (capping) of SRB stacking when an inadvertent ignition occurs. 
This event could happen in any of the four HB’s. Because of this, the results are superimposed on all HB locations 
to determine the final SSD. 



Figure 17: Sample case booster/plume configuration 


The SRB segments are assumed to be identical to the SSP segments for this sample case. It's assumed that the 
chamber condition of these incomplete stacks is approximately 100 psia and 5000° F: the chamber gas specific heat 
ratio was assumed 1.17. Nozzle geometry used in the SSP SRB's was assumed along with the same SRB/ET sitting 
position on the ML in the VAB. 

Because the boosters are uncapped, two plumes develop. The lower plume, because of the nozzle, has a larger exit 
Mach number (M c - 2) plume than the top plume (M c = 1). As a result, the lower plume sources are distributed 
farther out into the farfield and radiate sound further. In the near-field, the upper plume sound radiation is mainly 
confined to the VAB: little of its sound propagates to any significant area outside the VAB (see the figure below). 

OASPL Contour (dB) for VAB Nearfield 
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Figure 18: Near-field OASPL contour for sample case 




Upper Plutne Flow Direction (ceiling deflection) 
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In the farfield, the sound radiation is dominated mainly by the lower plume. After approximately 500 feet the 
OASPL decreases rapidly until it drops to insignificant levels past about 1000 feet (see the figure below). 



OASPL Contour (dB) for VAB Farfield (truncated to 0-100’ elevation) 



••c«0)dvv«tt IV) 

Figure 19: Farfield OASPL contour for sample case 

The SSD is drawn below with the SSP QD arc discussed earlier. The SSD in the sample case sits well within the 
previous SSP QD arc (see the figure below). 



Figure 20: VAB area SSD plot for sample case (OASPL lcc = 148 dB, OASPLopfo = 145 dB) 

From the transmission loss analysis, the LCC and OPF-3 impedances were not enough to reduce the OASPL below 
the 140 dB personnel limit; as a result these facilities would be unsafe as heavens. 

VII. Comments and Conclusions 

A methodology for determining a SSD for aeroacoustic hazards has been successfully developed. Future 
improvements of the method could aim to reduce the conservatism (likely in the range of 10-20 dB). A new module 
is planned such to model the effectiveness of ear protection equipment. Furthermore, a module to parametrically 
model the effectiveness of facility barrier shielding is likely. The addition of these two new' modules would add the 
dimension of mitigation to the model and methodology. 
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Introduction 

• In early 2007 efforts initiated to quantify inadvertent hazards in KSC 
processing areas for Constellation Program (CxP) 

• Use of the weight-based approach methodology resulted in very large arcs 

• Alternatively, hazards were laid out individually to be assessed 

• Some of the hazards considered were blast/fragmentation, radiative 
heating, aeroacoustics, toxics 

• End result was to determine a safe separation distance (SSD) for each 
individual hazard 

• The inputs, assumptions, criterion, and results of each individual hazard 
analysis were reviewed by peer committee 

• This presentation deals specifically with the Aeroacoustic hazard 
component 
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Assumptions 

• Burning stack/motor(s) are fixed in place and not pitched/tilted 

• Assumed VAB paneling and internal-dividing walls removed 

- simplifies analysis 

- provides most conservative acoustic fields for farfield 

• Plume(s) treated as axisymmetric 

- asymmetric spreading of the plume is not considered 

• Adjacent plumes in close proximity conjoined 






- only if plumes within 1-5 exit diameters 

- sound powers combined to create one plume 
Neglected atmospheric sound absorption 

- effect negligible for small regions considered 

- provides small degree of conservatism to analysis 



Neglected ground reflections (other than initial plume deflection) 
Assumed incident acoustic load on facility normal and symmetric 
Facilities simplified to rectangles for diffraction analysis 
Walls and windows treated as simply supported panels 
At frequencies where absorption data unavailable, assumed no absorption 



- difficult to find building material data below 125 Hz 

- majority of rocket noise sits in the 10 1 ~ 10 3 Hz area ^ 

- provides a large level of conservatism in analysis 
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Model Overview 



• Modified Eldred's Method 


Normalized Acoustic Power in Plume 

0* ' ' ? r r- 






Rocket plume divided into slices 
Each slice is treated a monopole source 



Receiver 




Flow Direction* 


Overall power of each slice dependant on 
location in plume 

Power maximum at ~1.5x core length 
(transitional region) 

Normalized converted to overall using, 

/Ax t (s)\ 

+ L w + lOlog ^ - J 


' \A .ps 


= 10 log 


X'Wtf 


W, 


OA 


Source Locations 
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Plume Source Directivity Index 


— s =0.001 


140 


160 


60 80 100 120 
Exhuast Relative Angle (deg.) 


180 


Each source has spectrum dependant on plume strength (in form on Strouhal number) 

Normalized converted to overall using, 

L f = 10 log (w(f,s)/W(s)p^^])+ L - 10 log (^SS±\ + lOlog (A f) 

' sa e ' ' sa amb ' 

Directivity of source with receiver dependant on Strouhal number and true angle with plume flow 
direction; SPL at frequency for given receiver and source, 

SPL p3 .f = ~ 10log(fi 2 ) - n + Di(p tff ) 

For any given receiver, sound pressure from each source summed at each frequency; resulting in SPL 

spectrum for receiver, ... 

r , fST SPL r*f 

SPL f = 10 log ( > 10 io 


Receiver SPL spectrum integrated to find OASPL 

OASPL = lOlog 


SPLf 
10 10 
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The diffraction-gain spectrum is computed from semi-empirical correlations using the building 

x = 108.354 m 


0.1 


L c /X - 


1 .0 

Minimum Dimension 
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OPF 3 



x -208 957 m 

x =254.139 m 


Woveiengfh 

The spectrums at the facility distances are then corrected using the diffraction-gain spectrum to find 
the acoustic excitation on the facility wall for the transmission analysis 


SPL'(f) = SPL(f) 

^ diffraction (/) 

Sample Acoustic Spectrum 
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• Based on the facility geometry and composition, wall 
transmission and internal absorption losses are calculated 

TL = 101og 10 (^) AL = 10 log 10 ) 

• Internal spectrum then computed by taking external diffracted 
spectrums and correcting for transmission/absorption losses 



materials 

125 HZ 

250 HZ 

500 HZ 

1000 HZ 

2000 HZ 

4000 HZ 

Painted Concrete: 

0 1 

005 

006 

007 

0.09 

008 

Carpet 

001 

005 

0 10 

020 

0.45 

065 

Window GUm: 

0 18 

0.06 

004 

003 

0.02 

0.02 

1 molrum Floor 

002 

003 

003 

003 

0.03 

0.02 

Ceiling Tilea: 

030 

051 

080 

0 85 

0.78 

066 

Gvpcum Board 

029 

0 10 

006 

0 05 

0.04 

004 

Remov able Wall* 

0 14 

0.10 

006 

004 

0.04 

0.03 


SPL' = SPV -TL + AL 

\ , i 


Diffracted Spectrum 




TLdiB) 



/. /, 0-V, /,, /, 
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Aeroacoustic Validation 



• Modified Eldred's 2 nd Method used for analysis 

• Original method (NASA SP-8072) has been updated with new core 
and directivity correlations 

• Super-position principle has been used to incorporate multiple non- 
conjoined plumes 

• Method used to predict acoustics around launch pad under STS 
launches 

• Method also used to predict launch acoustics for Ares-I/X 

• Sample data shown to right for Ares-IX test flight on October 28, 
2009 (data values not shown for ITAR/SBU considerations) 

• Excellent agreement between model and measurements in the near 
and midfield spectrums; ± 0.5 dB 

• Farfield spectrums from model over predict slightly in range of 
interest (10-1000 Hz); ± 7 dB 

• OASPL values (not shown for ITAR/SBU considerations) between 
model and measurements for nearfield, midfield, and farfield off by 
only + 1 dB 
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Personnel Exposure Limit 

• OSHA 1910.95(b)(2) would require a 115 dBA limit 

- OSHA table meant for industrial working environments and not inadvertent ignition events 

- Would result in overly conservative safe separation distances 

• Practical limit would be that associated with onset of physical impediment 

• Typical human ear threshold of pain occurs at 140 dB; ear drum rupture is typical at 160 
dB 

• Personnel exposure limit used was 140 dB (OASPL); taken at conservative head level (5') 



OSHA 1910.95(b)(2) 


Duration Per Day, in Hours 

Sound Level dBA, Slow Response 

8 

90 

6 

92 

4 

95 

3 

97 

2 

100 

1-1/2 

102 

1 

105 

1/2 

110 

1/4 or less 

115 

Note: When the daily noise exposure is composed of two or more periods of noise exposure of 
different levels, their combined effect should be considered, rather than the individual effect of each 
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Calculation Methodology 
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Aeroacoustic Procedure 
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Diffraction/Transmission Procedure 
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Sample Scenario 


• VAB is 4 th largest building in the world by volume 

• Consists of 4 High Bays (HB) used to process vehicles 
(Shuttle, Ares-I, Ares-V) 

• Ares-V vehicle uses (2) SRBs strapped to center stage 

• SRBs are assembled first in VAB on ML piece-by-piece 

• Sample scenario considers inadvertent ignition of (2) 
4.5-segment stacks in a single HB 

• Event could occur in any one of the HBs, so SSD 
calculated and placed at each HB center 

• Final arc drawn from centroid of HBs such that all 
SSDs encompassed 
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OPF-3 OASPL = 145 dB (UNSAFE) 



RED ARC = sample case GREEN ARC = SSP 


OASPL Contour (dB) for VAB Nearfield 


500 


Upper Plufne Flow Direction (ceiling deflection) 


1170 


1160 


400 


150 


o>300| 


(9 

C 


2001 


140 


h130 


1120 


100 1 


110 


Lower Pluhie Flow Direction (floor deflection) 


100 


100 200 300 400 

x-coordinate (ft) 


500 600 


F 


LCC OASPL = 148 dB (UNSAFE) 

OASPL Contour (dB) for VAB Farfield (truncated to 0-100’ elevation) 


i 

i 

i 
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Future Work 

• Obtain low frequency (< 125 Hz) data for facility materials 

- either by experimental testing or computationally (FEA) 

- remove some conservatism from analysis 

- could be as much as 2-3 dB for certain facilities 

• Develop ear plug/muff module using ANSI S12. 68-2007 to better determine mitigation 
levels 

- better determine effectiveness of safety protection 

- better quantify effectiveness of certain equipment for rocket acoustic spectrums (10-1000 Hz) 

• Parametric analysis of facility shielding 

- overlapping mitigation strategy from all hazard analyses involved blast shielding facilities in close 
proximity to VAB 

- analysis would determine effectiveness of shield size/construction to acoustic hazards 

• Develop GUI for aeroacoustics code and register in NASA database 

- easier to obtain funding to support work if registered 

- could make available to other agencies or possibly open-source (ITAR restricted) 
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BACKUP SLIDES 
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General Acoustics 
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SOUND POWER LEVEL 


(W\ 

L ~ = 10108 (wj 

where W is the acoustic power (typically in aero-acoustics some percentage of the mechanical power), 

and W 0 is the reference power (Watts) 

SOUND PRESSURE LEVEL 

SPL = 20 log (—) 

where p is the root mean square (rms) pressure, 
and p 0 is the reference rms pressure (Pa) 

p 0 is by standard 2x1 O' 5 Pa, this is equivalent to the approximate minimum 
pressure fluctuation detectable by the human ear 
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Source of sound 

Sound pressure 
(Pascal) 

Sound pressure level 
(dB re: 20 pPa) 

Theoretical limit for undistorted sound at 1 atmosphere 
environmental pressure 

101,325 Pa 

194.0937 dB 

Krakatoa explosion at 100 miles (160 km) in air 

20,000 Pa 

180 dB 

Simple open-ended thermoacoustic device 

12,000 Pa 

176 dB 

Ml Garand being fired at 1 m 

5,000 Pa 

168 dB 

Jet engine at 30 m 

630 Pa 

150 dB 

Rifle being fired at 1 m 

200 Pa 

140 dB 

Threshold of pain 

100 Pa 

130 dB 

Hearing damage (due to short-term exposure) 

20 Pa 

approx. 120 dB 

Jet at 100 m 

6 - 200 Pa 

110- 140 dB 

Jack hammer at 1 m 

2 Pa 

approx. 100 dB 

Hearing damage (due to long-term exposure) 

6*10 -1 Pa 

approx. 85 dB 

Major road at 1 0 m 

2xio- 1 -6x10' 1 Pa 

80 - 90 dB 

Passenger car at 1 0 m 

2x10-2- 2X10' 1 Pa 

60 - 80 dB 

TV (set at home level) at 1 m 

2x10- 2 Pa 

approx. 60 dB 

Normal talking at 1 m 

2xio-3-2xlO- 2 Pa 

40 - 60 dB 

Very calm room 

2xio- 4 -6x10- 4 Pa 

20 - 30 dB 

Leaves rustling, calm breathing 

6x10" 5 Pa 

10 dB 

Auditory threshold at 2 kHz 

2xl0 -5 Pa 
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A-Weighted Correction 



(/ 2 + 20.6 2 )y(/ 2 + i07.7 2 )(/~ 2 4- + 12200 2 ) 


.00 - 1 1 i 1 i 1 1 1 i 1 i. 1 1 1 1 j ! 1 1 — _ 1 1 1 1 1 1 

10 1 10 2 10 3 10 4 

f (Hz) 

- human ear more sensitive at certain frequencies 

- low frequency and high frequency sounds perceived to be not as loud as mid-frequency sounds 

- ear most sensitive to noise around 2-6 kHz range 

- A-weighted curve is the standard for quantifying sounds pressure levels 22 
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Aero-Acoustics 
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Decreasing Plume Deflection 
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Normalized Acoustic Power in Plume 
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increasing nozzle expansion ratio 



Dlifi.ff) = 10 log 


0.37 -+- 0.37 cas 4 (/7 ff ^) 


[(1 - 0.75cos(P' ff )y + 0.1687] ;5 [l + 310exp(-9/?,„ )] 


— 0.698log(s n (s)) — 1.67 


20 40 60 80 100 120 

Exhuast Relative Angle (deg ) 


140 


160 


180 


— S =0.001 

n 

— S =0.01 

n 

— S =0.1 
n 

— S =1.0 
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Plume Source Directivity Index 


l-i 
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Acoustic Diffraction 
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Diffraction Overview 

• Bending or spreading of waves around obstacle 

- larger wavelengths will bend around smaller objects or barriers, but not vice versa 

• Smaller wavelengths will not diffract around smaller barriers as effectively 

- result known as 'shadow zone' where the sound is quieter than elsewhere around object 

• Same conclusion can be made for long wavelengths/longer barriers, which result in 
a 'shadow zone' due to reflection 

• Result of diffraction is reflection, which can affect incident wave (enhance sound 
intensity levels on surface) 


Total Pressure = Incident Pressure + Pressure Reflected 
Ptot = Pj + 
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Average dB loss, [dB] 


Diffraction Analysis 

Wyle dB Loss Model 
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Minimum Dimension / Wavelength, [dimensionless] 




NOTE: Wyle model is average diffraction correction to be applied; if the value of L c / A is 

greater than 10, the dB gain is assumed to be 6 dB. 
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Acoustic Transmission 
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EFFECTS OF ACOUSTIC LOADS ON STRUCTURES 




Sound can either be reflected, absorbed, or transmitted through a wall 



ACOUSTIC ABSORPTION COEFFICIENTS 

MATERIALS 

125 HZ 

250 HZ 

500 HZ 

1000 HZ 

2000 HZ 

4000 HZ 

5” Thick Fiberglass 

0.50 

0.97 

0.99 

0.99 

0.97 

0.93 

Carpet on Concrete 

0.11 

0.14 

0.24 

0.40 

0.50 

0.60 

Heavy Window Glass 

0.18 

0.06 

0.04 

0.03 

0.02 

0.02 

Marble Tile 

0.01 

0.01 

0.01 

0.01 

0.02 

0.02 
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SOUND TRANSMISSION LOSS 
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• Transmission loss represents acoustic resistance of structure 

• Allows you to quantify the drop in sound through a wall 


SPL X % # SPL 2 

TL 


( 

SPL 2 = SPL X - 7X + 1 0 log 10 

V 


See J 


TL = - 101og IO (r) 


SINGLE WALL SOUND TRANSMISSION LOSS 


TL(dB) 



A \JA 
A A 


0.5 1 A 

A 


DOUBLE WALL SOUND TRANSMISSION LOSS 



Iog,o(/ >fc) 


log 10 (/) 
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